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Abstract — The  optical  form  birefringence  in  porous  silicon  films 
is  measured  by  analyzing  the  transmitted  interference  intensity  of 
a  polarization  interferometer.  A  novel  form  birefringence  model 
called  “boundary  condition  (BC)  model”  for  porous  materials  is  in¬ 
troduced  and  evaluated  experimentally  against  samples  of  porous 
silicon  films.  The  variation  of  optical  indexes  of  refraction  vs  the 
porosity  in  silicon  films  agrees  with  the  calculated  values  of  na/ ne 
within  1%  error  using  the  BC  model,  in  contrast  to  the  ~15% 
error  using  effective  medium  approximation  model. 

Index  Terms — Birefringence,  boundary  condition  (BC)  model, 
effective  medium  approximation  (EMA),  polarization  interferom¬ 
eter,  porous  silicon  (PSi),  simulation. 

POROUS  silicon  (PSi)  has  attracted  much  attention  and  the 
techniques  of  fabrication  and  characterization  of  PSi  have 
been  established  [1],  The  optical  birefringence  in  fabricated  PSi 
films  shows  discrepancy  between  the  calculations  using  existing 
models  and  the  experimentally  measured  ones  in  birefringence 
[2]-[4]  as  well  as  of  reflectivity  [5].  In  this  manuscript,  we  in¬ 
troduce  a  novel  “boundary  condition  (BC)  model”  and  evaluate 
the  performance  of  this  model  for  the  characterization  of  optical 
birefringence  in  PSi  films.  The  BC  model  not  only  accurately 
predicts  optical  indexes  ne  and  n0  from  the  structural  informa¬ 
tion  of  such  an  artificial  nanostructured  inhomogeneous  mate¬ 
rials,  but  also  traces  the  variation  of  optical  birefringence  vs  the 
porosity  (=  volume  of  air  pores/total  volume)  of  the  mate¬ 
rials,  making  it  useful  for  various  device  engineering. 

Due  to  the  cylindrical  symmetry  of  the  PSi  film  with  micro¬ 
pores  distributed  uniformly  on  the  surface  [4],  the  resulting  ar¬ 
tificial  form  birefringence  can  be  described  similar  to  that  of  a 
positive  uniaxial  anisotropic  crystal  leading  to  the  rotationally 
symmetric  index  ellipsoid  with  the  eigenvalues  (n0,  n0,  ne)  in 
the  principal  coordinate  system  (see  Fig.  1).  The  basic  analytic 
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Fig.  1.  Schematic  diagram  of  experimental  setup  and  the  index  ellipsoid  of  a 
Psi  film. 

equations  for  effective  permittivities,  which  is  analogous  to  the 
serial/parallel  capacitor  problem  in  a  circuit,  can  be  obtained 
easily  by  balancing  effective  energy  of  a  microscopic  volume 
including  just  one  air  pore  for  a  simple  BC  of  electric  field. 
Our  model,  for  the  sake  of  simple  derivation  of  analytic  index 
equations,  simplifies  the  shape  of  air  pores  by  assuming  square 
pillars  that  will  have  only  one  or  two  types  of  walls  leading  to 
BCs  parallel  and  perpendicular  to  the  electric  field  of  the  op¬ 
tical  wave.  Then,  £ejf  =  e^Ai/A)  +  £2(^2/^)  for  the  parallel 
BC  case,  and  l/£eff  =  (LijL)/ei  +  (Lo/L)/£2  for  perpendic¬ 
ular  BC  case.  A{=  A\  +  A2)  is  for  area,  and  L(=  L\  +  L2) 
for  the  length  [6].  The  value  of  extraordinary  index  of  refraction 
ne  for  PSi  can  be  determined  from  the  parallel  BC  approxima¬ 
tion  as  the  electric  field  is  parallel  to  z-axis,  and  thus  parallel 
to  all  boundary  walls  of  air  pores  (see  Fig.  1).  In  contrast,  the 
calculation  of  ordinary  index  na  requires  both  parallel  and  per¬ 
pendicular  types  of  walls,  as  the  electric  field  is  on  the  surface 
of  the  film.  The  analytic  equations  refined  for  form-birefringent 
dielectrics  are  given  as  follows: 

ne  =  'Jfi1  *  s' Air  +  (1  -  P2)  -esi  (2) 

where  (3  =  s/Ja. ir,  /air  is  the  porosity,  K  is  an  arbitrary 
correction  factor,  and  £a;r  and  £g;  are  the  permittivi¬ 
ties  of  air  and  silicon,  respectively.  The  correction  factor 
K  is  introduced  to  account  for  the  distortion  of  elec¬ 
tric  flux  without  significant  change  in  na,  and  K  will 
be  set  to  JFsimpieBC  =  1  for  a  “simple  BC  model,”  and 
K  =  Jvbc  =  /31/8  +  (1  —  Z?1/8)  •  2£Air/(£Si  ‘  [£Air  +  ^Si]) 
for  a  “BC  model.”  Notice  that  the  factor  Jvbc  approaches 
to  2£Air/(esi  •  [^Air  +  £si])  f°r  extremely  low  porosity 
as  /V16  ~  0,  otherwise  Kbc  is  also  nearly  1  because 
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/V16  ~  1.  Thus,  expanding  (1)  as  a  polynomial  of  /3 

leads  to  TIq  —  \Asi  /air  '  2  •  (£gi  ^-Air)/(^Si  T  ^Air) 
for  a  two-dimensional  (2-D)  dielectric  model  with  sparse 
cylindrical  air  impurity.  Note  that  a  typical  three-dimen¬ 
sional  (3-D)  dielectric  model  for  spherical  air  impurity  gives 
na  =  \Asi  -  /air  •  3  •  (eSi  -  £Air)/(2  •  £si  +  ^Air)  [6],  which 
is  valid  only  when  impurities  are  sparse  (/a„.  ~  0)  and  so  far 
apart  enough  not  as  to  affect  one  another.  In  the  following  we 
evaluate  the  validity  of  the  equations  for  ne  and  na  [i.e.,  (1) 
and  (2)]  in  the  case  of  dense  inhomogeneity  range  in  PSi  films 
with  various  porosity. 

PSi  films  are  fabricated  from  highly  doped  p-type  (100)  sil¬ 
icon  wafer  by  anodic  etching  in  aqueous  HF  (49  wt%)  solution 
with  ethanol  in  the  volume  ratio  of  3: 1  at  room  temperature  [7]. 
An  atomic-force  microscope  was  used  to  determine  the  average 
pore  size  of  less  than  20  nm  in  width,  whereas  the  porosity  was 
controlled  to  vary  for  different  samples  from  ~70%  to  ~85%, 
estimated  using  gravimetric  method.  As  explained  in  introduc¬ 
tion,  this  porous  film  has  optic  axis  collinear  with  the  direction 
of  the  pores.  Consequently,  the  effective  dielectric  constant  for 
optical  field  polarized  in  the  plane  of  incidence  (i.e,,  TM  or  ex¬ 
traordinary  polarization)  will  depend  strongly  on  the  angle  of 
incidence  [7].  A  schematic  diagram  of  our  experimental  setup 
for  characterization  is  shown  in  Fig.  1. 

The  wavelength  of  an  external-cavity  tunable  laser  source  of 
1.52  //m  is  chosen  below  the  band  edge  of  silicon  to  minimize 
optical  absorption  loss.  The  porous  silicon  sample  is  mounted 
on  a  precision  rotation  stage,  allowing  precise  control  of  the  sur¬ 
face  normal  orientation.  The  propagation  direction  of  the  optical 
beam  is  set  first  at  a  =  0  to  be  collinear  with  the  2  axis  (i.e,,  the 
pores)  of  the  PSi  film.  The  polarizer  and  analyzer  are  aligned 
perpendicular  with  respect  to  each  other  and  at  45°  with  the  x 
axis  in  the  x-y  plane.  When  the  angle  a  in  the  z-x  plane  is  in¬ 
creased  from  zero,  we  find  the  transmitted  intensity  at  the  output 
of  the  analyzer 

1(a)  =  ^  [a2  +  b2  —  2a.6cos(A(/)(Q;))]  (3) 

27r 

A 4>(a)  =  [ne(a)Le(a)  -  n0(a)L0(a )]  (4) 

where  a  and  b  are  the  amplitudes  of  transmitted  electric  fields 
for  ordinary  (perpendicular  to  the  plane  of  incidence,  TE)  and 
extraordinary  (in  the  plane  of  incidence,  TM)  polarizations,  re¬ 
spectively,  A 4>(a)  is  the  effective  phase  difference  for  ordinary 
and  extraordinary  polarized  components,  and  A,  ne^(a),  and 
Le(0)(a)  are  the  laser  wavelength  in  vacuum,  angle  dependent 
effective  refractive  indexes,  and  the  geometric  path  lengths  in 
a  film  for  extraordinary  (ordinary)  component  for  a  given  ro¬ 
tation  angle  a,  respectively.  Note  that  n0(a )  =  n0  (ordinary 

index )  for  all  a  and  ne(a)  is  1/ \J cos2  a/n%  +  sin2  a/n \  in  this 
setup,  where  ne  is  an  extraordinary  index  of  the  index  ellipsoid. 
The  polarization  interference  signals  are  obtained  by  measuring 
the  transmitted  optical  intensity  for  each  PSi  film  sample  [7]. 
This  intensity  variations  are  simulated  by  a  homemade  code  for 
this  birefringence  measurement.  The  measured  and  simulated 
polarization  interference  signals  versus  the  angle  of  incidence 
are  shown  in  Fig.  2. 
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Fig.  2.  Comparison  of  the  polarization  interference  signals  through  a  porous 
silicon  film  as  a  function  of  incident  angle  (a)  for  two  films  of  same  thickness 
of  228  fj, m  with  slightly  different  porosity  of  0.875  and  0.873  and  (b)  for  two 
films  of  similar  thicknesses  of  90  and  100  //  m,  but  with  quite  different  porosity 
of  0.86  and  0.76. 


Note  that  the  chirp-like  behavior  of  the  transmitted  intensity 
from  the  effective  phase  difference,  A helps  the  curve  fit 
process  to  determine  unique  optical  index  values  ne  and  n0  [7], 
which  are  varied  independently  to  find  the  best  intensity  curve 
fit.  The  index  pair  values  (ne,  nQ)  of  the  best  fit  and  the  corre¬ 
sponding  error  ranges  for  all  cases  are  summarized  for  compar¬ 
ison  in  Table  1.  The  intensity  curves  simulated  with  the  index 
pairs  at  upper/lower  values  in  Table  I  show  large  deviation  of 
out  of  tolerance  from  the  measured  ones.  Note  that  the  error 
range  decreases  as  the  film  thickness  increases  due  to  the  in¬ 
crease  of  phase  lag,  leading  to  more  phase  variation  data  for 
accurate  curve  fit. 

We  also  investigate  the  dependence  of  refractive  indexes  on 
the  film  porosity  for  various  models  [see  Fig.  3(a)].  The  re¬ 
fractive  indexes  from  using  the  BC  model  is  compared  with 
that  of  the  Bruggeman’s  effective  medium  index  approximation 
(EMA)  in  Fig.  3(b).  It  reveals  apparent  agreement  of  the  mea¬ 
sured  results  with  the  BC  model. 

The  inset  figures  in  Fig.  3(a)  represent  the  simplified  cross 
sections  for  microscopic  BCs  with  the  direction  of  electric  field 
denoted  by  an  arrow.  The  first  figure  from  the  left  is  for  770;iayered 
of  layered  structure,  the  second  one  for  770;sep.Si  of  isolated  Si 
column  structure,  the  third  for  na. bc  or  7t0,simpleBC  °f  isolated 
air  pore  structure  of  PSi  films,  and  the  last  one  corresponds  to  ne 
of  layered  or  pillar  structure  of  PSi  films  with  surface  parallel  to 
electric  field.  Note  that  n0tBC  is  similar  to  that  of  n0iSjmpieBC  f°r 
all  /a;r  (just  slightly  larger  than).  It  should  also  be  noticed  that 
the  curve  for  n0isep.Si  is  much  lower  than  that  for  n0iSimpieBC> 
although  the  only  difference  in  structure  is  the  connectivity  of 
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TABLE  I 

Comparisons  of  Optical  Indexes,  Error  Range,  and  Porosity  Values  for  Various  PSi  Film  Samples 


PSi  film 

Porosity 

(Gravimet.) 

Porosity 

(BCmodel) 

(fig,  Hq)  lower 

(ne,  Ho)best 

(ne,  n0)  upper 

(±D  nc,  ±0  n„) 

72  (gm) 

0.82 

0.85 

(1.5330,  1.3728) 

(1.6394,  1.4533) 

(1.7395,  1.5283) 

(±0.10,  ±0.08) 

90  (gm) 

0.71 

0.76 

(1.8641,  1.6248) 

(1.9235,  1.6668) 

(1.9812,  1.7079) 

(±0.06,  ±0.04) 

100  (gm) 

0.82 

0.86 

(1.5330,  1.3590) 

(1.6047,  1.4117) 

(1.6733,  1.4618) 

(±0.05,  ±0.05) 

198  (gm) 

0.82 

0.83 

(1.6733,  1.4635) 

(1.7066,  1.4877) 

(1.7393,  1.5117) 

(±0.033,  ±0.024) 

228  (gm) 

0.82 

0.873 

(1.5330,  1.3582) 

(1.5584,  1.3763) 

(1.5871,  1.3971) 

(±0.026,  ±0.019) 

228  (gm) 

0.82 

0.875 

(1.5330,  1.3541) 

(1.5515.  1.3673) 

(1.5692,  1.3803) 

(±0.018,  ±0.013) 

t  1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 
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Fig.  3.  Comparison  of  (a)  index  curves  from  various  models  as  a  function  of 
porosity  and  (b)  index  pair  coordinates  of  ( ne ,  nQ )  in  the  index  plane.  Four 
different  diagrams  in  (a)  denote  simplified  microscopic  geometric  arrangements 
of  pores.  The  arrow  represents  the  direction  of  the  electric  field. 

Si.  The  value  for  n0,i;«-eve<]  comes  from  an  unrealistic  BC  that 
all  boundaries  are  always  perpendicular  to  the  field,  only  obtain¬ 
able  for  layered  structure.  The  resultant  value  of  rrG, layered  curve 
is  the  lowest  for  any  /a;r  among  all  refractive  index  curves.  Fur¬ 
thermore,  sep.Si  is  larger  than  n0jayered  due  to  the  contribu¬ 
tion  of  the  increased  effective  Si  area  at  the  center  region  for  the 
same  /a;r,  as  shown  in  the  inset  figure  of  simplified  microscopic 
cross  section.  The  value  na  calculated  from  EMA,  which  satis¬ 


fies  n0  =  s/e^,  and  0  =  J2i  /,{(£;  -  e„r)/(c,-  +  2eav )},  devi¬ 
ates  much  from  the  measured  values,  and  it  is  larger  than  n0,BC 
for  low  porosity  (/a;r  <~  0.4)  and  smaller  than  n0j>c  for  high 
porosity  (/air  >  0.5) .  The  correlation  between  optical  indexes 

ne  and  na  for  various  porosities  can  be  viewed  more  easily  in 
Fig.  3(b).  It  is  apparent  that  the  BC  model  agrees  with  the  mea¬ 
sured  index  values,  while  the  curve  from  EMA  somewhat  devi¬ 
ates  from  the  measured  values.  The  introduced  BC  model  with 
simple  analytic  expression  provides  a  convenient  and  accurate 
prediction  for  the  optical  indexes  for  form  birefringence  in  PSi 
film. 

In  summary,  we  proposed  a  novel  optical  index  model,  called 
“BC  model,”  for  the  optical  form  birefringence  of  PSi  films, 
and  evaluated  it  experimentally  by  comparing  the  predicted  in¬ 
dexes  to  the  measured  ones.  The  optical  indexes  of  PSi  films  was 
characterized  accurately  using  optical  measurements  based  on 
the  data  from  polarization  interferometer.  The  optical  indexes  of 
PSi  films  are  predicted  much  more  accurately  by  the  BC  model 
than  any  other  pre-existing  models,  to  the  best  of  our  knowledge. 
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